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Helical magnetorotational instability of Taylor- Couette flows in the 
Rayleigh limit and for quasi-Kepler rotation 
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The magnetorotational instability (MRI) of differential rotation under the simultaneous presence of axial and azimuthal 
components of the (current-free) magnetic field is considered. For rotation with uniform specific angular momentum the 
MHD equations for axisymmetric perturbations are solved in a local short-wave approximation. All the solutions are 
overstable for B z ■ B^ 7^ with eigenfrequencies approaching the viscous frequency. For more flat rotation laws the 
results of the local approximation do not comply with the results of a global calculation of the MHD instability of Taylor- 
Couette flows between rotating cylinders. - With B^, and B z of the same order the traveling-mode solutions are also 
prefered for flat rotation laws such as the quasi-Kepler rotation. For magnetic Prandtl number Pm — * they scale with the 
Reynolds number of rotation rather than with the magnetic Reynolds number (as for standard MRI) so that they can easily 
be realized in MHD laboratory experiments. - Regarding the nonaxisymmetric modes one finds a remarkable influence 
of the ratio B^, / B z only for the extrema. For B$ B z and for not too small Pm the nonaxisymmetric modes dominate 
the traveling axisymmetric modes. For standard MRI with B z 3> B^, however, the critical Reynolds numbers of the 
nonaxisymmetric modes exceed the values for the axisymmetric modes by many orders so that they are never prefered. 
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1 Introduction 

It has been shown in previous publications starting with 
Hollerbach & Riidiger (2005) that the magnetorotational in- 
stability (MRI) under the presence of both current-free ax- 
ial and azimuthal components of the magnetic field ('He- 
lical' fields, hence HMRI) is always characterized by an 
eigenoscillation frequency. In combination with the vertical 
wavenumber the resulting instability pattern is thus an ax- 
isymmetric wave traveling along the rotation axis. In all our 
considerations the magnetic Prandtl number Pm plays the 
basic role. For Pm — » the HMRI scales with the Reynolds 
number Re rather than with the magnetic Reynolds number 
Rm as it does for the standard magnetorotational instabil- 
ity for an axial magnetic field. The questions arise whether 
this frequency reflects the geometry of the magnetic field 
and whether it is observable with real astrophysical ob- 
jects such as protoneutron stars and/or accretion disks. Also 
the relation to the Azimuthal MagnetoRotational Instability 
(AMRI, see Riidiger et al. 2007b) for (current-free) toroidal 
fields which is nonaxisymmetric and which scales with the 
magnetic Reynolds number for Pm — * must be consid- 
ered. We find the AMRI only weakly (if ever) influenced 
by the addition of an axial field which is not much stronger 
than the toroidal field. For Pm ~ 1 the difference between 
Re and Rm disappears so that the main differences of the 
instabilities also disappear. 



In the present paper we start with a local approximation 
using analytical methods for the most simple rotation law 
for constant specific angular momentum, i.e. the Rayleigh 
limit. Global calculations of the stability of the same rota- 
tion law between two rotating perfect-conducting cylinders 
and threaded by a helical current-free magnetic field lead to 
an almost perfect coincidence of the results of both meth- 
ods. This is no longer true, however, for more flat rotation 
laws such as quasikeplerian rotation in a finite gap where the 
differences of the short-wave approximation (which only 
holds for infinitely thin gaps) and global models are so 
strong that the results completely differ (Riidiger & Holler- 
bach 2007). 

2 Dispersion relation for very small gaps 

The dynamics of conducting fluids is described by the MHD 
equations 

du . . 1, 

(u ■ V)it = — V | p ■ 



at v 1 P 



EL 

2p 
——(B 



V)B + vA 2 u, (1) 
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0B 

— = rot (u x B) + 7jA 2 B, (2) 

and divtt = divB = 0. Here it is the fluid velocity, p the 
pressure, p = const the density and B is the magnetic field. 
In cylindrical symmetry this set of equations has the station- 
ary solution B = B c f,(R)e < p+B z e z where B z is a constant 
and B^ = I/R is current-free except on the rotation axis. 
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Following Lakhin & Velikhov (2007) we start with a 
local analysis. The perturbations are assumed to be axisym- 
metric and proportional to exp(7^ + \knR + ik z z), where 
7 = — ioj with (the real part of) ui as the eigenoscilla- 
tion frequency and k = (kn,0,k z ) the wave number vector. 
The local short-wave analysis is justified when k R ^> 1/R. 
Then the equations of small-amplitude perturbations of the 
steady-state flow u = fi(R)R are 
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Here u v = vk 2 is the viscous frequency, oj v = r\k 2 is the 
resistive frequency and the epicyclic frequency k is 



K 2 = 



1 d 



(9) 



P' is the total pressure perturbation. Then the dispersion 
relation 

k 2 

[(7 + u v )(i + w„) + ^|] 2 + -^k 2 x 

(7 + w^) +w A + i — - — (uj v -lu v ) 

k 2 
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results with 

(kB ) 2 



B 2 



POP 



J A<f, 



PopR 2 ' 



= (10) 



(11) 



The local dispersion relation is an equation of fourth or- 
der for the frequency. In the two limiting cases when either 
= or wa = the coefficients in Eq. (10) are real. In 
general, however, the coefficients of the dispersion relation 
(10) are complex. To find the stability criterion the method 
by Elstner, Rudiger & Tschape (1989) is applied. For ideal 
fluids Blokland et al. (2005) used a very similar approach. 



3 Rotation with constant specific angular 
momentum (the Rayleigh limit) 

For rotation with constant specific angular momentum we 
have k = which is called the Rayleigh limit in the 
Taylor-Couette community (see the MRI papers by Willis 



& Barenghi 2002; Rudiger, Schultz & Shalybkov 2003; Ve- 
likov et al. 2006). One obtains from Eq. (10) the critical 
angular velocity 

k 2 



F 



2k 2 \uja\ 

{uo\ + uj n uj v ) 2 + A(k 2 Jk 2 )uJ v uj v uo 2 A<j) 



1 1/2 



(w v + u v ) 2 u Arjl + {k 2 /k 2 )F(u% + uj, n Lu„) 



(12) 



with F = oj a + uj v uj v + 2(k 2 /k 2 )uJ A<f) . The frequency of 
marginally stable modes is 

LO = W A X 

(oj\ + uj n uj u ) 2 + 4(k 2 /k 2 )u> v w v GJ A(j> 



(ui v + u;„) 2 lj 2 a , + (k 2 /k 2 )F(uJ A + uj v oj u ) 



1/2 



(13) 



(Lakhin 2007, private communication). For a very weak az- 
imuthal magnetic field the instability limit is then defined 
by 

k 2 



n 2 = 



1 



4fc? 



(u A + uj ti uj v ) 2 + -j^^A<p [ 2 ( W A+ u v uj u ) - (uj v - LO v f 
Obviously, when the condition 

2(lV A + UJ^Uy) < (Ld v - UJ U ) 2 



(14) 



(15) 

is fulfilled the azimuthal magnetic field supports the insta- 
bility. Hence, the destabilizing effect of B^ disappears for 
Pm = 1 but it also exists for Pm > 1. 

With the free parameter a — k z /k and with 



aHa = Ha*, 



one finds from (12) and (13) 

2 w A0^A 



7Z = 2a 



F 



(16) 



(17) 



and F ~ lo v lo^ for weak magnetic fields. Then 
w = ^u^A <k L (lg) 

'' c LUv^l-j 

Hence, the oscillation frequency results from the simultane- 
ous existence of H, B z and B^ (Knobloch 1996). It exists 
despite of k = 0. The magnetorotational instability with 
helical but current-free magnetic fields is characterized by 
the existence of a frequency which only for the special case 
B<j, — vanishes (see also Blokland et al. 2005). As the si- 
multaneous existence of poloidal and toroidal components 
of the magnetic field is quite characteristic for cosmic ob- 
jects one could expect that the observation of the eigenfre- 
quency may serve as the proof of the existence of the MRI 
('magnetoseismology', see Blokland et al.). 
In the following the pitch angle 

(19) 
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(3/k z R and 



is used, and we write (3* 
Rc* = aRc = a — . 
Equations (12) and (13) then lead to 

Re ~ a, 



(20) 



(21) 



meaning that 

u> ~ uj v (22) 
for the marginal stable mode (Lakhin & Velikhov 2007). 
Obviously, even a small viscosity proves to be important for 
the system. Computations on the basis of the dispersion re- 
lation (25) under neglect of the viscosity (Liu et al. 2006) 
cannot lead to the same conclusion. With the renormaliza- 
tion 

= ~^=, (23) 
one finds 

w *~V^m, (24) 
independent of (3. At the Rayleigh limit it is thus the viscos- 
ity frequency which is emanated by the MHD system. 

We shall show that this relation indeed is valid for the ro- 
tation law fulfilling the Rayleigh condition, i.e. Q oc R~ 2 . It 
remains true for both local and global calculations. Though, 
for more flat rotation profiles (such as the Kepler law 
Q oc i?~ 3 / 2 ) the situation becomes more and more com- 
plicate, and the oscillation frequency starts to run with the 
Alfven frequency u A of the axial magnetic field (but not 
with the pitch angle (3*). Differences between the local 
and the global results become more and more large, only 
the Pm-dependence of the characteristic frequency remains 
weak. 



4 The numerical method 

A numerical method to solve the dispersion relation has 
been developed. The dispersion relation reads 

7 4 + ai 7 3 + a 2 7 2 + (a 3 + ib 3 )7 + «4 + i&4 = 0, (25) 
where 

ai = 2(lu v + uj v ), 

k 2 k 2 

a 2 = (lu v + lu„) 2 + 2(w| + lu v lu„) + -|k 2 + 4^|w| , 

k 2 

a 3 = 2(Ur, + UJ V )(U)% + LUjjLUv) + 2-^K 2 LU rl + 

k 2 

k 2 k 2 
a 4 = (lu 2 a + uj. n uj u ) 2 - A-^uj\Q 2 + -^k 2 {lo 2 a +lu 2 ) + 

k 2 

i>3 = — 8^|/2wawa<^, 

t „ k 2 2 k 2 z u a uj A4> 

h = -A-j^llLOA^A^-n + Uv) - k p — q — ^ _ UJ "'~ 

(26) 




10 -8 -|0" 6 10" 4 10" 2 10° 10 2 
Pm 

Fig. 1 Rayleigh limit: Critical Reynolds number vs. magnetic 
Prandtl number for axial fields (/3* = 0) and for helical fields 
(/3* > 0). The curves are marked with the oscillation frequencies 
cu* which obviously do not depend on the /3-values. 



The dispersion relation (25) is now solved numerically for 
k = 0. To this end the (complex) frequency is again nor- 
malized with ^lu„oj v . Then from (26) 




a 2 = (VP^+ -jLJ) +2(1 + Ha* 2 ) +4/T 2 Ha* 2 , 
V VPm / 

a 3 = 2 (1 + Ha* 2 ) ( VPm+ 




+ 4/3* 2 Ha* 2 ( VPm + 



04 = (1 + Ha* 2 ) 2 - 4Ha* 2 Re* 2 Pm + 4/3* 2 Ha* 2 , 
b 3 = -8/3*Ha* 2 Re*Vplii, 
6 4 = -4/3*Ha* 2 Rc*(l + Pm). 



(27) 



For given Pm, (3* and Ha* the dispersion relation can be 
solved in the following way. For too small Re* all the so- 
lutions of (25) will have negative ^(7) so that the pertur- 
bations decay. For a certain Re*, however, the first real 
part of one of the roots becomes positive so that the dis- 
turbance grows and the flow becomes unstable. Then that 
Ha* is searched which provides the lowest Re*. This min- 
imum Re* is shown in Fig. 1 for various /?* in its depen- 
dence on the magnetic Prandtl number. Also the resulting 
normalized frequencies uo* are given which indeed do not 
depend on (3* . Note the basic difference of the solution for 
(3* = and (3* > 0. Clearly for Pm ^ the solution 
with (3* = scales with Rm* = Pm • Rc* while the so- 
lutions with [3* > scale with Re. For Pm = 1 no basic 
difference exists between the solutions as there is no differ- 
ence between Re and Rm. The consequence of the different 
scalings for small magnetic Prandtl number is a rather small 
critical Re* for (3 > 0. 

One can also take from Fig. 1 that for (3* = Re oc 
l/VPm. We know this relation as a particularity of the 
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Fig. 2 Rayleigh limit: Critical Reynolds numbers vs. toroidal 
field component /3* for fixed magnetic Prandtl number (Pm = 
1CT 5 ). On the vertical axis: Re* = 317, Ha* = 1. 
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Fig. 3 Rayleigh limit: Hartmann number and oscillation fre- 
quency for fixed magnetic Prandtl number (Pm = 10~ J ) for very 
strong toroidal field components. Note the constancy of ui* over 
the whole range of Reynolds numbers and pitch angles f3* > 0. 



Rayleigh limit 1 which changes to Re oc 1 /Pm (i.e. Rm ~ 
const) beyond the Rayleigh limit (see Rudiger & Hollerbach 
2004). 

We now ask for the function Re* = Rc* (/?) for fixed 
and small Pm. Figure 2 shows the characteristic profile for 
Pm = 10 -5 . A minimum Reynolds number exists for (3* ~ 
0.5 which is smaller by two orders of magnitude than Re* 

for /3* =0. 

When the Reynolds number decreases for growing (3 the 
frequency uj* increases to the value \/Pm and remains then 
constant (see Eq. 24). For small Pm it does not depend on (3 
and/or Re*. The amplitude confirms the basic result (22). 

Beyond the minimum Ha* oc 1//3* and Re* oc (3* 
which can be observed in Fig. 2. 

In summary, the MRI for the Rayleigh limit J? oc R~ 2 
exists for not too high ratios of azimuthal and axial mag- 
netic field components. It is basically oscillating but not for 
the exceptional case of B,p = 0. The oscillation frequency 
approaches the viscous frequency uj v . Pitch angels of order 
ten strongly destabilize the MRI, there is a deep and wide 



1 x10 6 



52 5x10 5 




It follows directly from the relation 124 = for fi = re = 0. 



Fig. 4 Kepler flow: The characteristic minimum in the Reynolds 
number profile disappears. The toroidal field component stabilizes 
the MRI. On the vertical axis: Re* = 66670, Ha* = 259. Pm = 

icr 5 . 



trough in the profile Re = Re(/3) as shown in Fig. 2. Note 
that for j3 of order ten the scaling with Re for Pm — > still 
exists which is characteristic for HMRI. 

5 Kepler flow: local solutions 

In order to model the flat rotation law of a Kepler flow we 
have to work with (26) and k = Q. The numerical results 
are shown in Fig. 4. They strongly differ from the results for 
the Rayleigh limit. The characteristic trough in the Reynolds 
number profile in Fig. 2 disappears. The toroidal field al- 
ways stabilizes the MRI. After Fig. 4 the Reynolds number 
grows with /?*, i.e. 

Re* ~ 2 • 10 5 /3*. (28) 

There is no instability anymore with very small magnetic 
Reynolds number or small Lundquist number. Obviously, in 
the local approximation the solution in Fig. 1 which scales 
for Pm — > with Re instead of Rm disappears for too flat 
rotation profiles. 

For the resulting frequency uj* Fig. 5 provides uj ~ 
0.8wa independent of the pitch angle (3. In units of the ro- 
tation frequency, however, 

7? K W- (29) 

This result depends on the pitch angle of the magnetic field 
but does not depend on the magnetic Prandtl number. 

6 Global TC containers 

With a local approximation we have considered the instabil- 
ity of differential rotation under the presence of a current- 
free magnetic field with helical geometry, i.e. B Z B^ 7^ 0. 
The local approximation only describes flows with an ex- 
tremely small radial extension. The instability for such field 
geometries forms an overstable pattern of axisymmetric per- 
turbations of flow and field. The corresponding eigenoscil- 
lation only vanishes for B$ — 0. Beyond the Rayleigh limit 
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20 



Fig. 5 Kepler flow: The Hartmann number Ha* and the oscil- 
lation frequency to* are almost constant for /3* > 1 and they are 
almost equal so that again u> ~ 0.8wa- Pm = 10~ 5 . 



Lakhin & Velikhov (2007) find in the short-wave approxi- 
mation that the eigenfrequency fulfills ut ~ an without any 
relation to the magnetic field. 

The basically unknown quantity in local approximations 
is the resulting wave number a. We shall switch, therefore, 
to the consideration of global models where the wave num- 
ber no longer is an unknown quantity. We shall find that the 
eigenoscillation u of HMRI equals the viscosity frequency 
also for global solutions for all pitch angels (3 but only at 
the Rayleigh line. For more flat rotation laws the resulting 
eigenfrequency of the HMRI modes proves to be fixed by 
the Alfven frequency of the vertical field. 

For comparison with the results of the local approxima- 
tion we first consider a container with a small gap between 
the cylinders (i?i n /-R ut = 0.67) which are assumed as per- 
fect conductors. The rotation law between the cylinders is 



with 

_ m -fj 2 n 



l-fj 2 



where 



n 



HQ 



b = 



p2 Q 

-"-in "in j 



(30) 



(31) 



(32) 



Rout Qi. 

f2i n and tl out are the imposed rotation rates of the inner 
and outer cylinders with their radii R m and i? ut- The 
Rayleigh limit is reached if fin = if , i- e - = 0.44 
for the considered small-gap container. The axial field B z 
is uniform while the toroidal field is current-free in the 
gap, i.e. B^/Bz = (3R 1TL / R.The. Hartmann number and 
the Reynolds number are the dimensionless numbers of the 
problem, 

B z Rq 



Ha 



Re 



i?; n i?Q 



(33) 



where R = (Ri n (Rout — Rin)) 1 ^ 2 is taken as the unit of 
length. We have used R^ 1 as the unit of the wave number 
and H in as the unit of frequencies. We shall also use the 
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Fig. 6 Reynolds number (top) and eigenfrequency (normalized 
with J?in, bottom) for a small-gap container (fj = 0.67) filled with 
gallium (Pm = 10~ 6 ) vs. the rotation rate ratio /in- The vertical 
axis at the left represents the Rayleigh limit. The curves are marked 
with their (3- values. Note the clear dependence of uj/{2 in on the 
pitch angle (5. 



magnetic Reynolds number Rm = Pm • Re. Details of the 
numerics can be found in Riidiger et al. (2007a,b) and ref- 
erences therein. 

In Fig. 6 the Reynolds number (top) and the normalized 
oscillation frequency u> / 17; n (bottom) are given for very 
small magnetic Prandtl number (Pm = 10~ 6 ) and for vari- 
ous values of (3. The main result is that for larger j3 the crit- 
ical Reynolds decreases and the normalized (!) oscillation 
frequency grows. In consequence, the product of Reynolds 
number and oscillation frequency ui / J?; n should not vary to 
much with the /3-values (see Fig. 7, top). At the Rayleigh 
limit we indeed find one and the same value for the consid- 
ered product independent of (3 leading back to (21) yielding 
the basic relation (22). At the Rayleigh line the oscillation 
frequency is thus exactly given by the viscosity frequency. 
This is not true, however, beyond the Rayleigh line. There 
we have the phenomenon that the oscillation frequency nor- 
malized with the viscosity frequency strongly depends on 
the /3-value (Fig. 7, top). However, the /3-dependence be- 
yond the Rayleigh line is reduced if the Hartmann number 
of the axial field is used (Fig. 7, bottom). For (medium) 
(3 ~ 4-8 the curves show Re • (cj/O in ) oc Ha independent 
of p. This means 



(34) 
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Fig. 7 The same as in Fig. 6 but for the product of Reynolds 
number and eigenfrequency (top) and this product divided by Ha 
(bottom). Note that at the Rayleigh line the product of Reynolds 
number and eigenfrequency does not depend on f3 and note that 
outside the Rayleigh line for medium (3 the use of the magnetic 
normalization is successful. 

almost independent of [3 for medium values of (3. A weak 
dependence on jift remains. Equation (34) means that for 
given magnetic Prandtl number the amplitude of the vertical 
field determines the oscillation frequency. This relation is 
not valid for too small nor too large values of j3. For too 
small (3 the influence of the toroidal field is too weak and 
for too large (3 the influence of the axial field is too weak. 

7 Kepler flow: global solutions 

In contrast to the previous results we have demonstrated 
with nonlocal calculations that even for Kepler rotation 
HMRI occurs if at least one of the boundaries is sufficiently 
conducting (Rudiger & Hollerbach 2007). The short-wave 
method used in the foregoing sections does not reflect these 
results. In the following the Taylor-Couette calculations for 
flat rotation laws have been repeated i) to probe the basic re- 
sults with another code and ii) to find more details about the 
dependence of this effect on the magnetic Prandtl number 
and about the eigenfrequencies and iii) to include nonax- 
isymmetric modes into the considerations. In order to ob- 
tain exact solutions we remain in the frame of the Taylor- 



10-8 10-6 10" 4 10-2 10° 10 2 
Pm 

Fig. 8 The critical Reynolds numbers for quasi-Kepler rotation. 
Note that the fields with (3 = and f3 = 10 for Pm -> scale 
with Rm and the fields with f3 = 2 and f3 = 4 scale with Re. 
Perfect conducting cylinders, m = 0, fj = 0.5. 

Couette flows so that only quasi-Kepler rotation laws can 
be considered. A quasi-Kepler rotation law fulfills at both 
cylinders the Kepler law Q oc R~ 15 so that the condition 
fiQ = fj 15 results. Here (as also in Rudiger & Hollerbach 
2007) with the standard gap of fj = 0.5 is worked so that 
/tin = 0.35 holds for the quasi-Kepler approximation. 

The experiment PROMISE works with a rotation law 
with /l^ = 0.27. The experimental setup utilizes the exis- 
tence of the trough in the relation of Re and (3 where the crit- 
ical Reynolds number of the instability becomes very low. 
For more flat rotation laws one expects an increase of the 
critical Reynolds number. Lakhin & Velikhov (2007) have 
shown that for n = fl an instability exists, with frequencies 
about one order of magnitude smaller than the basic rota- 
tion, but the trough for medium /3-values disappears. All the 
solutions are now scaling with Rm rather than with Re. We 
shall show, however, that this result does not reflect the sit- 
uation in Kepler disks which obviously cannot be described 
in a local approximation. 

Figure 8 gives the critical Reynolds numbers as a func- 
tion of (3 and Pm. It is clearly shown that for small (3((3 = 0) 
and large ((3 = 10) the instability scales with Rm while 
for Pm — > the instability for (3 = 2 and (3 = 4 scales with 
Re. This is the usual constellation for HMRI. For Kepler 
rotation there is no basic difference close to the Rayleigh 
line. For Pm = 1 no distinction exists between Re and Rm 
so that all differences disappear between the solutions for 
different (3. The oscillation frequencies u> / Q ln are of order 
0.1 (with a weak /3-dependence for small Pm) and they only 
slightly depend on the magnetic Prandtl number (see Fig. 9, 
top). One finds 

fV 4 ,l (35) 
J7j n Ha 

for Pm between 10~ 6 and 1 (Fig. 9, bottom). Hence 

Pm 1/4 , (36) 



© 2008 WILEY- VCH Verlag GmbH & Co. KGaA, Weinhdm 



www.an-journal.org 



Astron. Nachr. / AN (2008) 



7 



3 




10 -8 10"6 10" 4 10" 2 10° 10 2 

Pm 



Quasi-Kepler 



en 

LlJ 
QQ 



< 

X 




10-a 10" 6 10" 4 10- 2 10° 10 2 
Pm 





1 0.00 



1 .00 



0.10 



0.01 




1 0° 



Fig. 9 The normalized eigenfrequency (u)/Oin, top), the critical 
Hartmann numbers (middle) and the numerical basis of Eq. (35) 
for quasi-Kepler rotation (bottom). In the bottom plot the curves 
are for /? = 2, 4, 6. Parameters as in Fig. 8. 



independent of (i. Only the amplitude of the axial field fixes 
the eigenfrequency. The influence of the magnetic Prandtl 
number is small. This is a numerical result for not too small 
Pm, the physical meaning of the Pm-factor is not yet com- 
pletely clear. The result (36) for Pm = 10~ 5 is smaller by 
one order of magnitude than the result derived from the dis- 
persion relation. 

7.1 Nonaxisymmetric modes 

Also nonaxisymmetric solutions with cxp(im</>) have been 
computed with this model which represent solutions drifting 
along the azimuth. Riidiger et al. (2005) find the solution 
of this problem for m = 1 for various fiQ but only for a 



Fig. 10 The critical Reynolds numbers for standard MRI (i.e. 
/? = 0) vs. magnetic Prandtl number at the Rayleigh limit (top) and 
for quasi-Kepler rotation (bottom). The curves are marked with 
their mode number m. In both cases for small Pm the differences 
are huge. 



fixed magnetic Prandtl number. Always the mode with m = 
1 needs a higher Reynolds number to be excited than the 
mode with m — 0. At the Rayleigh limit and for Pm = 
10~ 5 the critical Reynolds number proves to be 3-10 6 . 

7.1.1 Standard MRI 

For standard MRI (i.e. (3 = 0, see Riidiger & Zhang 2001; 
Ji, Goodman & Kageyama 2001) the calculations have 
been repeated for various magnetic Prandtl numbers and for 
m = 0,1, 10 (Fig. 10). For all nonvanishing m the critical 
Reynolds numbers lie above the curve for m = 0. It is also 
understandable that this effect grows for smaller Prandtl 
numbers. The differential rotation increases the magnetic 
dissipation of nonaxisymmetric patterns growing with m 2 ?? 
and T) runs with 1 /Pm. This situation exists for both steep 
rotation law (top) and also for flat rotation laws such as 
the Keplerian one (bottom). In general, the nonaxisymmet- 
ric modes are strongly stabilized, by the smoothing action 
of differential rotation but also by the appearance of the 
azimuthal drift as argued by Mikhailowskii et al. (2008). 
Standard MRI is obviously a basically axisymmetric phe- 
nomenon. 
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Fig. 11 The same as in Fig. 8 but with the nonaxisymmetric 
mode m — 1 included. The solid lines give the axisymmetric 
modes (m = 0) for various while the dashed line belongs to 
the nonaxisymmetric (m = 1) modes for all 0. 
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Fig. 12 The same as in Fig. 1 1 but for the Hartmann number. 



7.1.2 HMRI 

The motivation for the inclusion of nonaxisymmetric modes 
for helical magnetic geometry is that for large the toroidal 
field dominates the axial field so that previous computations 
concerning the m = 1 instability for azimuthal (current- 
free) magnetic fields with oc 1/R ( AMRI) should be 
concerned. In Fig. 1 1 the dashed line gives the instability 
limit for = 10. One finds that it scales with Rm for Pm — > 
also known for AMRI (see Rudiger et al. 2007b). 

For the mode with m — 1 the critical Reynolds num- 
ber for small Pm exceeds the critical Reynolds number for 
tyl = so that the axisymmetric traveling-mode is much 
easier to excite than the nonaxisymmetric solution. The op- 
posite is true for Pm > 10~ 3 . In this case, for > 10, the 
nonaxisymmetric modes are more easy to excite. 

We also computed the nonaxisymmetric solutions for 
= 2 and = 4. We find that their instability limits can 
not be distinguished from that for = 10. It is clear that 
also the solutions with > 10 comply with the curve for 
= 10. The Hartmann number, however, becomes smaller 
and smaller for growing which, however, is not a surprise 
as for azimuthal MRI (also for the Tayler instability) the 
toroidal field strength is important, i.e. the product Ha. 



Note that for small Pm the Hartmann numbers for m = 
1 are much higher than for m = (Fig. 12). This may be 
a simple consequence of the fact that AMRI scales with the 
Lundquist number S = \/Pm Ha. 

7.2 Results 

The basic results from our calculations for quasi-Keplerian 
rotation laws are: 

1 . For < < 4 and Pm < 1 the axisymmetric HMRI has 
always the lowest Reynolds number. The traveling-wave 
solution with m = which is observed in the PROMISE 
experiment forms the instability which is most easiest to 
excite also for Kepler rotation. Opposite statements (Liu 
et al. 2006) cannot be confirmed. 

2. The oscillation frequency for > equals the viscos- 
ity frequency but only at the Rayleigh line. For quasi- 
Kepler rotation it proves to be the Alfven frequency of 
the vertical magnetic field. 

3. For larger (0 ~ 10) and Pm < 10~ 3 the lowest 
Reynolds numbers belong to the axisymmetric HMRI 
while for Pm > 10~ 3 they always belong to the nonax- 
isymmetric AMRI. 

4. The neutral instability line for m = 1 does not depend 
on 0, i.e. the AMRI is not concerned by the existence of 
axial fields. 

The calculations demonstrate the strong influence of the 
magnetic Prandtl number for the MHD instability theory. 
The formula SMRI (axisymmetric, stationary, scaling with 
Rm) + AMRI (nonaxisymmetric, drifting, scaling with Rm) 
= HMRI (axisymmetric, oscillating, scaling with Re) exists 
only for small Pm. 
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